Energy loss of charged particles moving parallel to a magnesium surface 
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We present ab initio calculations of the electronic energy loss of charged particles moving outside 
a magnesium surface, from a realistic description of the one-electron band structure and a full 
treatment of the dynamical electronic response of valence electrons. Our results indicate that the 
finite width of the plasmon resonance, which is mainly due to the presence of band-structure effects, 
strongly modifies the asymptotic behaviour of the energy loss at large distances from the surface. 
This effect is relevant for the understanding of the interaction between charged particles and the 
internal surface of microcapillaries. 

PACS numbers: 71.10.Ca, 71.45.Gm, 73.20.At, 78.47.+p 



I. INTRODUCTION 

Collective excitations at metal surfaces (surface plas- 
mons, SP) 1 - are well known to play a key role in a wide 
spectrum of science, ranging from physics and materials 
science to biology^ Here we focus on one specific situa- 
tion where surface plasmons play a key role: the energy 
loss of charged particles moving near a metal surface. 
This work has been partially motivated by recent the- 
oretical and experimental studies of the interaction of 
highly charged ions^ with the internal surface of micro- 
capillaries and nanocapillaries^^ whose interpretation 
calls for an understanding of the asymptotic behaviour 
of the energy loss at large distances from the surface. 
Our ab initio calculations of the electronic energy loss 
of charged particles moving outside a magnesium surface 
indicate that the finite width of the SP resonance, which 
is mainly due to the presence of band-structure effects^ 
strongly modifies this asymptotic behaviour. 

When a charged particle is placed in front of a metal 
surface, the distribution of electrons in the surface and 
the bulk is modified, an induced charged density is built 
up, and both single-particle and collective excitations 
are created which arc ultimately responsible for the elec- 
tronic energy loss of the external charged particle. Mov- 
ing charged particles can also lose energy as a result of 
the interaction with the nuclei of the solid; however, this 
energy loss is negligible compared to the electronic energy 
loss, unless the probe charge moves at velocities that are 
extremely small compared to the Fermi velocity vp of the 
solid. 

For a weak interaction between the external charged 
particle and the electrons of the metal surface, the elec- 
tronic response can be treated within linear-response the- 
ory. In the case of charged particles moving inside a solid, 
nonlinear effects arc known to be crucial in the inter- 



pretation of energy-loss measurements however, non- 
linear corrections have been shown to be less important 
when the charged particle moves outside the solid.— On 
the other hand, the electronic response of the metal sur- 
face is expected to be strongly dependent upon the actual 
structure of the ground-state electron density. Therefore, 
an accurate description of the electronic energy loss of 
charged particles moving near a metal surface lies mainly 
in the understanding of two basic ingredients: the elec- 
tronic properties of the ground state of the target and 
the linear response of a many-electron system to exter- 
nal perturbations. 

Existing self-consistent calculations of the energy-loss 
spectra at solid surfaces are based on either the jcllium 
model 1 - 1 - or a one-dimensional (ID) model potential that 
still assumes translational invariance along the surface.— 
Here we report the first three-dimensional (3D) ab ini- 
tio calculation of the electronic energy loss that incorpo- 
rates the full band structure of the solid surface. Specif- 
ically, we calculate from first principles the energy loss 
of charged particles that move parallel to the (0001) sur- 
face of Mg. We use the random-phase approximation 
(RPA)^ and we focus on the limit in which the moving 
particle travels at large distances (compared to the Fermi 
wavelength) from the surface. We demonstrate that even 
for a free-electron-like metal such as Mg band-structure 
effects on the finite width of the SP resonance strongly 
modify this asymptotic behaviour. 

The paper is organized as follows. First of all, in 
Sec. II we describe our full ab initio treatment of the 
wavevector- and frequency-dependent electronic response 
of valence electrons, based on a realistic description of the 
one-electron band structure. From the knowledge of the 
density-response function we derive the surface-response 
function of the system, whose imaginary part is related 
to the rate at which an external potential generates elec- 
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tronic excitations at the solid surface. In Sec. Ill, we de- 
rive an explicit expression for the electronic energy loss 
of charged particles moving parallel to a solid surface 
(stopping power of the solid), which we obtain from the 
knowledge of the imaginary part of the GW self-energy 
of many-body theory.— In Sec. IV, we present numeri- 
cal calculations of the surface-response function and the 
electronic stopping power of a Mg(0001) surface. We 
compare our results with the stopping power obtained 
with the use of either the jellium model or a ID model 
potential. In Sec. V, our conclusions are presented. Un- 
less otherwise is stated, we use atomic units throughout, 



i.e., e 



1. 



where A is a normalization area, q is a 2D wave-vector 
in the surface Brillouin zone (SBZ), and g and g' denote 
2D reciprocal lattice vectors. For an external potential 
of the form 



A 



(4) 



the rate w(u) takes the forrn^ 

SBZ 



(5) 



II. SURFACE-RESPONSE FUNCTION 

Take a system of N interacting electrons of density 
no(r) exposed to a frequency-dependent external poten- 
tial V ext (r; uj). Keeping terms of first order in the exter- 
nal perturbation and neglecting retardation effects, the 
rate w(u>) at which the frequency-dependent external po- 
tential V ext (r;uj) generates electronic excitations in the 
many-electron system^ is found to be given by the fol- 
lowing expression: 

w(uj) = -2Im J dx J dr'V ext (r; uj) x(r, r'; uj) V ext (r'; uj), 

(1) 

where x(r, r'; uj) represents the so-called density-response 
function of the many-electron system. In the 
framework of time-dependent density-functional theory 
(TDDFT)fi£ the exact density- response function of an 
interacting many-electron system is found to obey the 
following Dyson-type equation: 



X(r,r';w) = x°(r,r';o>) + / dr t / dr 2 X° (r, r t ; uj) 

x {u(n,r 2 ) + r c K](r 1 ,r 2 ;w)}x(r 2 ,r';w). (2) 

Here, x°(r, r'; u>) is the density-response function of non- 
interacting Kohn-Sham electrons, u(ri,r2) is the bare 
Coulomb interaction, and / IC [rio](r,r';w) represents the 
so-called exchange-correlation (xc) kernel, which is the 
second functional derivative of the xc energy functional 
evaluated at the ground-state electron density no(r). In 
the RPA, f xc [n ](r,r';uj) = 0. 

In the case of a periodic surface, we introduce 2D 
Fourier expansions of the form 



SBZ 



where w(q, uj) denotes the rate at which the external po- 
tential generates electronic excitations of frequency cu and 
2D wave vector q: 



Air 

w{ci, uj) = — Imff(q,cj), 



(6) 



with 



ff(qiw) = J dz J dz' e q{z+z>) Xg=o :g '=o(^,z';q,w). 

(7) 

In the RPA, the Fourier coefficients \ Stg i (z , z' ; q, u>) are 
found to obey the following matrix equation: 

Xg.g' [z,z'\ q, uj) = Xg, g ' {z,z'\ q, uj) 

+ ^2 / dzi / dz 2Xg, gl (^, ^i;q,w) 
gl J 

x v gl {z 1 ,z 2 ;a L )x^ 1 ^{z2,z'- 1 (\,uj), (8) 

where v g (z, z';q) denote the 2D Fourier coefficients of 
the bare Coulomb interaction v(r, r'), 



v s (z,z';q) 



2% „-|q+gl|z-*'l 



q + g 



(9) 



x(t,t',uj) = I ]T E e4(q+s) ' r " e " (q+g ' ) ' r ' l x g , g '(^^ / ;q^), 

(3) where 



and Xg g ' i z i z '\ q^ uj) represent the Fourier coefficients of 
the density-response function of noninteracting Kohn- 
Sham electrons. 

For positive frequencies, the imaginary part of the 
Fourier coefficients x g g '( z ? z '; q ; ^-0 is easily obtained 
from the spectral function S g g i(z, z'\ q, uj), as follows 

lmxl !S '(z,z';q,uj) = -irS gtg ,(z, z';q,u), (10) 



q s,s' 
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S°, g , q,u,) = 7 EE E (^,n(r)|e- i ^ r ii|^ +q ,„Kr))(V'k +q ,nKrO|e i ( q+8 ') r »|^,„(r')) 
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X <5(ek,n - £k+q,n' + w)- (H) 



Here, the sum over n and n' run over all occupied and 
unoccupied bands, respectively, and £k,n and ipk,n(z) rep- 
resent, respectively, the single-particle energies and wave 
functions of a Kohn-Sham Hamiltonian with an effective 
potential that is periodic in the plane of the surface. For 
the evaluation of the real part of Xg g '( z i z \ Q; lJJ )^ we P er ~ 
form a Hilbert transform of the corresponding imaginary 
part,iZ. 

The function <?(q, oj) of Eq. ([7]) is the so-called surface- 
response function, which represents a key quantity in 
the description of both single-particle and collective ex- 
citations at solid surfaces, and whose imaginary part 
[see Eq. ((BJ] yields the rate at which an external poten- 
tial generates electronic excitations. Equation {7J shows 
that only the diagonal Fourier coefficient Xg,g( z ) z '\ 1) w ) 
with g = enters the expression for the surface- 
response function. Nevertheless, the full matrix na- 
ture of Xg, g '( z : z ''i 1i w ) an d Xg,g' i z > z '\ <L w ) is involved 
when solving Eq. ([8]). These arc the so-called local-field 
effects^ which typically play an important role in the 
presence of strong electron-density inhomogeneities but 
are found to be negligible in the case of simple metals 
like Mg^ 



III. ELECTRONIC STOPPING POWER 

Let us consider a probe particle of charge Z\ moving 
in an inhomogeneous many-electron system. The decay 
rate t" 1 of the particle in the state <f>i(r) with energy 
is obtained from the knowledge of the imaginary part of 
the self-energy S(r,r';£ ?; ), according toiS 

rr 1 = -2 [ dr [ dr>*(r)Im£(r, r'; ^(r'). (12) 



In the GW approximation of many-body theory^ and 
replacing the probe-particle Green function by that of a 
noninteracting particle, one finds: 



1 = -1Zl \ dr J dv'4>*{v) 



x r f (r')lmW(r, r'; £l - £/ )^(r)0 4 (r'), (13) 
/ 

where the sum is extended over a complete set of final 
states 4>f( r ) ofenergye/. W(r, r'; £j — £/) is the screened 
interaction of the system, which is related to the inter- 



acting density-response function as follows 

W(r,r';u) = v(r,r') 

+ I dvi I dr 2 v(r, ri)x(ri,r 2 ; w)w(r 2 ,r'). (14) 



In the case of a recoilless point particle moving at a 
given impact vector b with nonrelativistic velocity v, the 
probe-particle initial and final states can be described by 
plane waves in the direction of motion and a Dirac 5 
function in the transverse direction, i.e., 



1 



e iV V<$(r±-b), 



(15) 



where rj_ represents the position vector perpendicular to 
the projectile velocity. One then finds that the decay rate 
of Eq. (|13|) can be written as follows 



(16) 



where T is a normalization time and P q is given by the 
following expression: 



dq' 



,ib-(q+q') 



47T 2 f°° 

x imVF(q, q', w)S(w - q ■ v)6(u> + q' • v). (17) 

Here, f2 is a normalization volume and W(q, q';w) rep- 
resents the double Fourier transform of the screened in- 
teraction W(r, r';ui). 

The quantity P q entering Eq. (fT6|) can be interpreted as 
the probability for the probe particle to transfer momen- 
tum q to the many-electron system. Hence, the stopping 
power of the many-electron system, i.e, the average en- 
ergy lost by the particle per unit path length is found to 
be given by the following expression: 



dE 1 ^ 

- = z £(q.v)P q 



(18) 



where L is a normalization length and q- v represents the 
energy transfer associated to the momentum transfer q. 

Now we restrict our attention to the case of charged 
recoilless particles moving with constant velocity v along 
a definite trajectory at a fixed distance z from a periodic 
solid surface. If one introduces 2D Fourier expansions of 
the form of Eq. ©, then Eqs. $T7D and JTHJ) yield the 
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following expression for the stopping power: 
dE 



dx 



07 2 SBZ 
iZj \ JK b 



vA ^ ^ 

g.K q 



q- v 



x ImW g) g +K (2,z;q,q- v), (19) 



the sum over K being restricted to those reciprocal lat- 
tice vectors that are perpendicular to the velocity of the 
projectile, i.e., K • v = 0. 

At this point, we focus on the special situation where 
the coordinate z is located far from the surface into the 
vacuum. Equation Q14p shows that under such condi- 
tions the Fourier coefficients W SyS ' (q, u>) take the follow- 
ing form: 

Wg, g '(z,2;q,w) =v g (z,z,q}S g!S > 
2nq 



q + g q- 



- ffg , g .(q,c) e -^+sl+|q+ g 'l), ; (20) 



where 



2tt 

ff g ,g'(q,w) = — - 



I q+g' I *' 



(21) 



which for g = g' = yields the surface-response function 
entering Eq. ([6]). 

The symmetry of the one-particle Bloch states results 
in the following identity: 



(22) 



with S representing a point group symmetry operation 
in the periodic crystal. As a consequence, the stopping 
power of Eq. (fl"9|) can be evaluated from the knowledge 
of the screened interaction corresponding to wave vectors 
lying in the irreducible element of the surface Brillouin 
zone (ISBZ). If crystal local-field effects are neglected 
altogether, introducing Eq. (|2"Tj) into Eq. (fT!)]) yields 

q s 

^-(5q-v)e- 2 ' Sq ' z Im 3 (5q,q.v), 

where <?(q, a;) represents the surface-response function of 
Eq. ®. 

In the simplest possible model of a solid surface, in 
which a semi-infinite electron gas described by a Drude 
dielectric function e(ui) = 1 — u 2 /w 2 is separated by a 
planar interface at z = from a semi-infinite vacuum, 
both Eqs. (fT!?|) and l|2"5|) reduce, for particle trajectories 
outside the solid (z > 0), to the classical expression^ 



■). (23) 



(z) = Z 1 2 ^-X o (2w s zA0, 



(24) 



where Kq(x) in the zero-order modified Besscl function^ 
and u> s is the SP frequency uj s = u) p /\/2. This expression, 
which is known to hold for high particle velocities (v >> 
vf) at large distances from the surface (z » Xf) shows 
that under these conditions the energy loss is dominated 
by the excitation of surface plasmons. 



IV. RESULTS 

Magnesium (ls 2 2s 2 2p 6 3s 2 ) is a monoatomic solid with 
the hexagonal close-packed crystal structure. The in- 
put of our parameter-free first-principles stopping-power 
calculations is the surface-response function </(q, uf) of 
Eq. ([7]), which we have calculated for the (0001) surface 
of Mg. The results presented below have been found to be 
well converged for all velocities under study. The single- 
particle Kohn-Sham orbitals ^k,n(r) entering Eq. (fTTj) 
were expanded in a plane-wave basis set with a kinetic- 
energy cutoff of 13 Ry. In the Fourier expansion of the 
noninteracting density-response matrix \ s g ' (<7z > 9 Z 1 Qi w ) 
we included 101 g z -vectors and the components with 
g = g' = only because of the negligible local-field ef- 
fects along the surfaced In Eq. (fTTj) all occupied and 
unoccupied energy bands up to 50 eV above the Fermi 
level were taken into account. Numerically, in the eval- 
uation of S® g /(.<7z, g' z , q, u>) the 5-function was replaced 
by a Gaussian with a broadening parameter a = 0.1 cV. 
The sampling of the SBZ required for the evaluation of 
both the surface-response function of Eq. (J7J) and the 
stopping power of Eq. (|23[) has been performed including 
7812 mesh points in the SBZ. We set Z\ = ±1, but our 
results can be used for arbitrary values of Z\, as the stop- 
ping power is, within linear-response theory, proportional 
to Z 2 . 

We compare our first-principles calculations with the 
results that we have also obtained by replacing the ac- 
tual (0001) surface of Mg by (i) a jcllium surface with an 
electron-density parameter r s — 2.66 (corresponding to 
an electron density equal to that of valence electrons in 
Mg) and (ii) a model surface described by the ID poten- 
tial of Ref . [22|. This potential describes the main features 
of the surface electronic structure of Mg: the energy gap 
and the Shocklcy surface state at the f point (k = 0) 
of the SBZ; in this case, we have used films of up to 41 
layers of atoms with a lattice parameter a — 4.923 a. u. 
corresponding to a film thickness of 100.92 a. u., and the 
work function has been taken to be $ = 3.66 eV. 

Our first-principles calculations employ a supercell ge- 
ometry with slabs containing 16 atomic layers of Mg that 
are separated by vacuum intervals. The slab geometry 
imposes a low limit for the momentum q below which 
the SP splits into two slab excitations of the form=. 



to±(q) 



,(lTe-« L ) 



qL-sl/2 



(25) 



with L representing here the slab thickness. This draw- 
back can be softened by increasing the slab thickness. 
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Figure Q] shows the self-consistent calculations of 
the imaginary part of the surface-response function, 
Img(q, u>), that we have obtained from Eq. ((7]) in the 
RPA for (i) a semi- infinite jellium surface (dashed lines), 
(ii) the ID model surface potential of Ref. [22| (thin solid 
lines), and (iii) the actual (0001) surface of Mg (thick 
solid lines). For the low 2D wave vectors q under study, 
the energy-loss spectra are clearly dominated by a SP 
contribution at w s ~ 8eV, which seems to first shift to 
lower frequencies, as q increases, and then, from q ~ 0.02 
on, towards higher frequencies. This figure shows that for 
the small values of q considered here both jellium and ID 



model calculations (dashed and thin solid curves) over- 
estimate the SP energy; our calculations show, however, 
that for larger values of the 2D wave vector these simpli- 
fied models predict accurate values of the SP energy. 

The important message of Fig. [1] is that at small val- 
ues of q the actual linewidth of the SP is considerably 
larger than that obtained with the use of ID jellium- 
like models. This important effect is expected to impact 
considerably the actual stopping power of the solid sur- 
face, especially at large distances from the surface where 
the energy loss is dominated by the excitation of surface 
plasmons associated to very low wave vectors. 




co (eV) co (eV) 

Figure 1: Imaginary part of the RPA surface-response function of Mg(0001), Ini(/(q|| , uj), as a function of u>, for various values 
of the magnitude q of the wave vector. The thick solid line represents self-consistent first-principles calculations, the dashed 
line represents the corresponding results obtained for a semi-infinite jelli um surface, and the thin solid line represents the 
corresponding results obtained by using the ID model potential of Ref. |23 . 



Now we focus on the special situation where an exter- 
nal recoilless particle of charge Z\ = 1 is moving with 
constant velocity along a definite trajectory at a fixed 
distance z far from the surface into the vacuum. Fig. [2] 
exhibits the stopping power of Mg(0001) for this moving 
particle, as obtained from Eq. (|2"5)) in the RPA for (i) 
a semi- infinite jellium surface (dashed lines), (ii) the ID 



model surface potential of Ref. [22| (thin solid lines), and 
(iii) the actual (0001) surface of Mg (thick solid lines). At 
the velocities under consideration (v > vp), the energy- 
loss spectrum of charged particles moving far from the 
surface into the vacuum is dominated by long-wavelength 
surface excitations (small q), i.e., by the excitation of sur- 
face plasmons. Hence, we might be tempted to conclude 
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z (a.u.) z (a.u.) 

Figure 2: Stopping power of Mg(0001), versus z, for a recoilless particle of charge Z\ = 1 moving with various velocities. The 
thick solid line represents self-consistent first-principles calculations. The dashed and thin solid lines represent the corresponding 
results obtained for a semi-infinite jellium surface and by using the ID model-potential of Ref. I22I respectively. The dotted 
lines represent the classical energy loss of Eq. (f24)l . 



that Eq. (represented in Fig. [2] by a dotted line) 

should be a good representation of the actual stopping 
at z >> Xp. Indeed, Fig. [5] shows that this classical 
limit is in excellent agreement at large values of z with 
the results obtained with the use of ID jellium-like mod- 
els (dashed and thin solid lines). However, it is important 
to note that these models do not account for the intrinsic 
linewidth of surface plasmons which, as a result of inter- 
band transitions that arc absent in these simplified mod- 
els, dominates the energy loss at large distances. As the 
velocity increases (see the lower panels of Fig. [5]), lower 
values of the wave vector enter the excitation spectrum 
leading to an increased impact of the intrinsic surface- 
plasmon linewidth on the stopping power and, therefore, 
more pronounced differences between the stopping power 
of a jellium-like surface (dashed and thin solid lines) and 
the real surface (thick solid lines). 

Figure [3] shows the ratio between the first-principles 
stopping power of the real Mg(0001) surface and the stop- 
ping power of the corresponding jellium surface [which 
at large distances coincides with the classical result of 
Eq. ([H])] for the four values of the velocity considered 
in Fig. [21 The ab initio energy loss is considerably larger 
than that obtained from the classical Eq. ([24]) at all large 
distances from the metal surface. As z — ► 00, the clas- 
sical Eq. (J2U) (which assumes the linewidth of the sur- 
face plasmon to be zero) decays exponentially; indeed, 
at these large distances from the surface the stopping 
power is dominated by the finite intrinsic linewidth of 
the surface plasmon, leading to a ratio that increases ex- 
ponentially. This exponential increase also occurs when 



N 




Figure 3: (Color online) Ratio 

{-dE/dx) abmltto {z)/(-dE/dxY emurn (z), versus z, for a 
recoilless particle of charge Zi = 1 moving with various 
velocities outside theMg(OOOl) surface. In the inset, larger 
distances z from the surface are considered. 



the stopping power of the real surface is replaced by that 
of a semi-infinite electron gas described by a Drude di- 
electric function of the form e(u>) = 1 — uc/u}(us + 17), 7 
being a damping parameter that would give account ap- 
proximately for the finite surface-plasmon linewidth^ 
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V. SUMMARY AND CONCLUSIONS 

We have carried out first-principles self-consistent cal- 
culations of the surface-response matrix and the stopping 
power of the (0001) surface of Mg. Our results indi- 
cate that band-structure effects (and. in particular, in- 
tcrband transitions) play a key role in the description of 
the asymptotic behaviour of the stopping power far from 
the surface, even in the case of a free-electron-like metal 
such us Mg. In particular, we find that the lincwidth of 
the surface-response function is considerably enhanced 
at small wave vectors, which yields an increased energy 
loss of charged particles moving far from the surface that 
cannot be described by simplified ID jcllium-likc models. 
This important effect, which should be present in the case 



of all metal surfaces, is expected to be relevant for the 
understanding of the interaction between charged par- 
ticles and the internal surface of microcapillaries. New 
experiments along these lines would be desirable. 
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